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Density engineering of an oscillating soliton/ vortex ring in a Bose-Einstein condensate 

Itay Shomroni, Elias Lahoud, Shahar Levy, and Jeff SteinhaueiB 

Technion-Israel Institute of Technology 

When two Bose-Einstein condensates (BEC's) collide with high collisional energy, the celebrated 
matter wave interference pattern results. For lower collisional energies the repulsive interaction 
energy becomes significant, and the interference pattern evolves into an array of grey solitons. The 
lowest collisional energy, producing a single pair of solitons, has not been probed. We use density 
engineering on the healing length scale to produce such a pair of solitons. These solitons then 
evolve periodically between vortex rings and solitons, which we image in-situ on the healing length 
scale. The stable, periodic evolution is in sharp contrast to the behavior of previous experiments, in 
which the solitons decay irreversibly into vortex rings via the snake instability. The evolution can 
be understood in terms of conservation of mass and energy in a narrow condensate. The periodic 
oscillation between two qualitatively different forms seems to be a rare phenomenon in nature. 

PACS numbers: 03.75.Lm, 03.75.Kk, 47.32.cf, 47.37.+q 
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We consider two condensates separated in space, and 
given an initial collisional energy E co \ , which could be in 
the form of kinetic energy, or potential energy in a mag- 
netic trap. The condensates are then allowed to collide in 
the z direction. For negligible interactions, the two col- 
liding condensates will interfere like two non-interacting 
plane waves, forming a standing wave interference pat- 
tern with wavefunction ip oc cos(fcz), where the wave- 
length of the pattern is A = 2ir/k : and k is given by 
h 2 k 2 /2m — E co \/N where m is the atomic mass and N is 
the number of atoms. The density of the pattern is there- 
fore of the form p], [2] n = \ip\ 2 oc cos 2 (kz). The initial 
energy E co \ supplies the kinetic energy of the standing 
wave, proportional to d 2 il)/dz 2 , as well as the negligible 
interaction energy Nji/A. However, if E co \ is decreased to 
become comparable to iV/x/4, the form of the fringes will 
be modified by interactions. It will become energetically 
favorable to form an array of grey solitons, rather than 
a cosine-squared pattern [2|]. The criterion E co \ < iV/i/4 
can be written as A > 47r£, where £ is the healing length, 
the approximate width of a grey soliton [3|, [4| . In other 
words, as E co \ decreases, the wavelength of the standing 
wave interference pattern increases until the wavelength 
reaches the size of a soliton. At this point, the fringes will 
evolve into solitons. A cosine-squared fringe and a soliton 
are actually very similar in form, as seen in Fig. [T^l, b. 
As E co \ decreases further, the number of solitons will 
decrease, maintaining conservation of energy [5]. The 
lowest collisional energies will result in two grey solitons 
traveling in opposite directions. Such a low collisional 
energy is achieved in our experiment by initially sepa- 
rating the two condensates by a barrier whose width is 
on the order of £, as shown in Fig. The energy of 
this initial configuration relative to the ground state is 
on the order of E co \ = AT/i(£/L), where L is the length of 
the condensate in the z direction. This is comparable to 
the energy of a grey soliton computed in one dimension, 
which is on the order of 

E sol = N^/L) (1 - v 2 /c 2 f 2 , (1) 
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FIG. 1: Solitons and matter- wave interference fringes, 
(a) The order parameter near the transition from interference 
fringes to solitons. The bold solid curve depicts an interfer- 
ence fringe with A = 47r£. The narrow solid curve shows a 
soliton with v = 0. Both of these order parameters are real. 
The dashed (dash-dotted) curve indicates the real (imaginary) 
part of a grey soliton, with v = 0.53c. (b) The density as- 
sociated with the curves of part (a), (c) The initial density 
and phase of our experiment, created by density engineering, 
(d) The density and phase of two counterpropagating solitons. 



where v is the subsonic soliton velocity, and c is the speed 
of sound in the condensate [3]. E co \ is thus sufficient to 
create a pair of solitons. 

Solitons can also be created by phase engineering @, 0] . 
There is a phase step < 7r across the grey soliton, 
given by Acj) = 2cos~ 1 (v/c). By imposing this phase 
step, the density profile will evolve into that of a soli- 
ton. We can understand this process by considering the 
current density J = (h/m)nV<j>. By imposing the phase 
step, this current density is significant in the region of the 
step only. This corresponds to a divergence in the flow 
field, V J = dJ z /dz ^ 0. Since the condensate obeys 
the continuity equation V J = —dn/dt, The divergence 
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depletes the density, forming the soliton minimum. 

We use density engineering to form solitons, a com- 
plementary technique to phase engineering. We impose 
the density profile, and the phase evolves into that of a 
soliton. The process is illustrated in Fig. [T]l which shows 
two solitons separated by a distance d. They are moving 
in opposite directions, with the corresponding opposite 
phase steps. Our initial configuration shown in Fig. [Tfc 
can be thought of as the d = case of Fig. [T]i. Our ini- 
tial configuration thus corresponds to two overlapping, 
counterpropagating solitons. 

In a one-dimensional (ID) condensate, solitons are sta- 
ble objects. In three dimensions however, grey solitons 
were seen to decay irreversibl y in to vortex rings via the 



snake instability 
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which was 
We find a 



originally measured in optical solitons [1 
very different, stable behavior shown in Fig. [2] by a sim- 
ulation of the 3D Gross-Pitaevskii equation (GPE) 
The soliton evolves into a vortex ring, which subsequently 
evolves back into a soliton. This process repeats itself 
many times with a period of 7-8 ms. We call this a "pe- 
riodic soliton/ vortex ring". The periodic soliton/ vortex 
ring is stable, surviving the reflection from the end of 
the condensate, as well as the subsequent collision in the 
middle of the condensate. The periodic soliton/ vortex 
ring has a constant phase step (0.77T for our experiment), 
as shown in Fig. [2], k. 

First, we will understand the evolution of the soliton 
into a vortex ring in terms of conservation of energy. Due 
to the higher density near the center of the condensate, 
the soliton moves fastest in this region [7]. This causes 
the soliton to curve 0, as shown in Fig. [2F. Due to the 
curvature, the area of the soliton has increased, corre- 
sponding to an increase in the soliton energy. To con- 
serve the energy of the soliton, the increase in area must 
be compensated by an increase in v, by equation [H This 
corresponds to a reduction in the depth of the soliton, 
as shown in Fig. [TJd. The soliton thus moves with ever- 
increasing speed and decreasing depth, until the density 
minimum vanishes altogether, leaving a vortex ring, as 
shown in Fig. Eb, g, and i. The appearance of the vortex 
ring does not violate Kelvin's theorem (conservation of 
circulation), because the vortex core can enter and exit 
the condensate via the nearby boundary [l7j |. 

The fate of vortex rin gs i n a quantum fluid is a long- 



standing field of study [17|, llS[. In superfluid He at 



non-zero temperatures, it was suggested that the radius 
of a vortex ring shrinks in time due to dissipation, until it 
evolves into a roton. We can understand the very differ- 
ent evolution of our vortex ring into a soliton in terms of 
conservation of mass. We consider a vortex ring in an in- 
finite homogeneous condensate 0, [l9[ , as illustrated [l8[ 
in Fig. [21. The density is constant in time, maintained 
by the zero divergence of the flow (V • J = 0). In other 
words, after the flow traverses the center of the vortex 
ring from right to left, it returns to the right outside the 




FIG. 2: Simulation of the periodic soliton/vortex ring. This 
simulation shows 3 x 10 4 atoms, (a)-(i) Slices of the conden- 
sate. The shades of grey indicate the density, (a) The conden- 
sate during the turn-off of the barrier, (b), (e), (h) The soliton 
stage of the periodic soliton/vortex ring, (c), (g), (i) The vor- 
tex ring stage. The axis of the 2 counterpropagating vortex 
rings is horizontal in the image. The slice of the vortex ring 
appears as two dark spots, separated vertically, (d), (f) The 
transitions between the soliton and vortex ring stages, (j) The 
shades of grey indicate the phase, and the red vectors indicate 
J. The image is an enlargement of the area between the red 
lines of part (c). The values and — 0.7tv indicate the phase 
on either side of the phase step, (k) As in part (j), except that 
the enlargement refers to part (e). (1) An illustration of the 
slice of a vortex ring in an infinite, homogeneous condensate, 
(m) An illustration of the 2 counterpropagating vortex rings 
of part (c), showing a surface of constant density (blue) as 
well as a surface of constant flow speed (red). 



vortex ring. In contrast, our vortex ring is close to the 
edge of the condensate, as shown in Fig. [2k, g, and i. 
The flow lines therefore have no return path to the right 
side, causing V - J ^ 0. The divergence depletes the den- 
sity, forming the soliton minimum. This situation is very 
similar to that of phase engineering of a condensate de- 
scribed above. In fact, the vortex ring has a phase step, 
as shown in Fig. [2J. The phase step occurs on a length 
scale Az on the order of the radius of the vortex ring, 
which is also the radius R of the condensate. Thus, for 
creating solitons from vortex rings, R should not be too 
much larger than £. This can also be understood energet- 
ically. The periodic soliton/vortex ring requires that the 
soliton and the vortex ring have the same energy. The 
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FIG. 3: In-situ images of the vortex rings and solitons. The 
experimental parts of the figure show averages of 2 to 6 de- 
structive phase contrast images, employing small detuning. 
The condensates contain 1 x 10 4 atoms, (a) The initial con- 
dition with the barrier on. (b) The height of the barrier as a 
function of time (not to scale). The solid curve corresponds 
to density engineering of solitons. The dashed curve corre- 
sponds to the creation of sound pulses, (c) Simulation of the 
initially formed pair of solitons at r = 2 ms. This is an in- 
tegrated image which includes finite resolution and depth-of- 
field. (d) In-situ image of the two count erpropagating solitons 
at 2 ms. (e) Simulation of the vortex ring stage at r = 5 ms. 
This is an integrated image which includes finite resolution 
and depth-of-field. (f) In-situ image of two counterpropagat- 
ing vortex rings at 5 ms. (g) Simulation of the soliton stage 
at r = 9 ms. This is an integrated image which includes fi- 
nite resolution and depth-of-field. (h) In-situ image of two 
counterpropagating solitons at r = 9 ms. 

ratio of these energies contains a factor of placing 
an upper limit on R. By the GPE simulation, we find 
that the vortex ring evolves into a soliton for R < 13£. 
This requirement implies that the number of atoms in 
the condensate should not be too large. For our trap 
frequencies for example, simulations show the periodic 
soliton/ vortex ring for N < 7 x 10 4 . 

Our experimental apparatus is basically described in 
Ref. [20]. After reaching BEC, the trap frequencies are 
adiabatically reduced by half in 350 ms to 117 Hz and 
13 Hz in the radial and axial directions, respectively. The 
resulting condensate contains 1 x 10 4 or 3 x 10 4 87 Rb 
atoms, as specified in the captions of Figs. [3] and [H The 
condensate with 3 x 10 4 atoms has \ijh — 580 Hz (h is 
Planck's constant), corresponding to a minimum healing 
length of £ = 0.32 /im, which is found at the center of 
the condensate. The healing length is larger at the lower 
densities found away from the center. The radii of this 
condensate are 89£ and R = 10£ in the axial and ra- 
dial directions respectively. The condensate with 1 x 10 4 
atoms has even larger £ and smaller R. 

Density engineering is performed by a laser beam, blue 
detuned by 6 nm from the 780 nm resonance. The con- 
densate is split in the radial direction by the highly elon- 
gated laser beam with an axial diameter of 4.4£ (1.4 //m). 
By the ID GPE simulations of Ref. [5], as well as our 3D 
GPE simulations, this diameter of roughly 4£ is the max- 
imum which will create a single pair of solitons. 
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FIG. 4: The speed of the periodic soliton/vortex ring. The 
condensates contain 3 x 10 4 atoms, (a) Integrated profiles of 
the periodic soliton/vortex ring. Each profile is an average 
of several images, except for the ms profile. The dashed 
curves are the result of the simulation, (b) Closed circles in- 
dicate the position of the periodic soliton/vortex ring as a 
function of time. Open circles indicate the position of the 
sound pulses. The solid and dotted curves are simulations 
of the periodic soliton/vortex ring and sound pulses, respec- 
tively, (c) The sound pulses. Each profile is an average of 
several images. The dashed curves indicate the simulation, 
(d) Two counterpropagating sound pulses at 5 ms. 



We adiabatically ramp up the laser barrier in r ramp = 
200 ms to a height U = ft(3000Hz), creating a BEC with 
a density minimum, as shown in Fig. [3^. The potential 
is then rapidly turned off in 1 ms, resulting in the initial 
condition illustrated in Fig. This sequence is indi- 
cated by the solid curve in Fig. [3)3. As seen in Fig. [3]i, 
one can distinguish the two counterpropagating solitons 
after 2 ms, which result from the interference between 
the two condensates of Fig. [3^, colliding with low en- 
ergy. Figure [3]i corresponds to Fig. [T]i. Figure [3t shows 
the simulation at 2 ms with imaging effects, including 
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integration of the density perpendicular to the image, 
finite resolution, and finite depth of field. The experi- 
ment (Fig. [3]i) is seen to agree well with the simulation 
(Fig-Efc)- 

After a time r = 5 ms, we obtain in-situ images of the 
pair of counterpropagating vortex rings of the periodic 
soliton/ vortex ring, as shown in Fig. [3j. In the image, 
the two vortex rings are seen as weak density minima. 
Again, good agreement is seen with the simulation of 
Fig. [3^. After each of the two vortex rings has evolved 
back into a soliton, we image the soliton stage at r = 9 ms 
as shown in Fig. [3fr. Two strong, slit-like solitons are 
seen. We find that the solitons are more visible for the 
small atom number shown in the figure. Good agreement 
with the simulation of Fig. [3^ is seen. 

In order to measure the speed of the periodic soli- 
ton/vortex ring, we observe its position as a function 
of time, as shown in the integrated profiles of Fig. [4^, 
and the filled circles of Fig. [4)3. We see that the speed 
(the slope) is approximately constant in time, and is 
therefore the same for the soliton stage and the vortex 
ring stage. The speed is observed to be v = 0.66 ± 
0.08 mms -1 , for N = 3 x 10 4 . The simulation gives 
the same results qualitatively and quantitatively, as in- 
dicated by the dashed and solid curves of Fig. 0^, b. 
This speed is of the same order of magnitude as the 
speed v ~ (h/2mR) ln(1.59i?/^) « 0.3 mms -1 predicted 
for a vortex ring in an infinite, homogeneous conden- 
sate SEEll. 

The speed of a grey soliton is less than c, so it is 
interesting to compare the speed of the periodic soli- 
ton/vortex ring to c. We measure c by the technique 
of Ref. [2l|. Specifically, we rapidly turn on the bar- 
rier in 100 fis and leave it on as indicated by the dashed 
curve of Fig. [3)3, creating two counterpropagating sound 
pulses, as shown after 5 ms in Fig. [Hi. The profiles of 
Fig. [4k and the open circles of Fig. [4)3 show the posi- 
tion of the sound pulses as a function of time, giving 
c — 1.24 ± 0.07 mms -1 , which agrees with the theoreti- 
cal value [22] of c = vV 2m = 1.15 ± 0.09 mms -1 . The 
dotted and dashed curves of Fig. [4)3, c show the results 
of the GPE simulation, which agree well with the experi- 
ment. We thus find that the periodic soliton/ vortex ring 
moves slower than the speed of sound, at v = 0.53c. 

In conclusion, we have achieved density engineering of 
a pair of counterpropagating grey solitons, which are the 
result of a low-energy collision between 2 BEC's. The 
close relationship between solitons and matter wave inter- 
ference fringes elucidates the quantum mechanical nature 
of solitons in a BEC. Due to the inhomogeneous nature 
of the narrow condensate, each soliton evolves into a pe- 
riodic soliton/ vortex ring. We explain this evolution in 
terms of conservation of mass and energy. Since vortex 
rings and solitons can be considered to be quasiparti- 
cles d, Qj], perhaps the periodic soliton/ vortex ring 
could additionally be described as a Rabi oscillation be- 



tween these quasiparticle states. Such a Rabi oscillation 
is seen in optics, between two soliton states [23|. As a 
macroscopic non-linear entity, the periodic soliton/ vortex 
ring has the unusual property that it oscillates between 
two qualitatively different forms. This presents a puz- 
zle to all branches of science, to categorize this effect by 
finding other systems with an oscillating nature. 
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